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Abstract 

We introduce an anomalous top quark coupling (right-handed current) into Standard 
Model Lagrangian. Based on this, a more complete calculation of b — ► S7 decay including 
leading log QCD corrections from ra top to Myy in addition to corrections from My/ to 
mi, is given. The inclusive decay rate is found to be suppressed comparing with the case 
without QCD running from mt to M\y except at the time of small values of e.g. 
when /Jj = —0.08, it is only 1/10 of the value given before. As |/f | goes smaller, this 
contribution is an enhancement like standard model case. From the newly experiment of 
CLEO Collaboration, strict restrictions to parameters of this top-bottom quark coupling 
are found. 
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1 Introduction 



The standard model(SM) has achieved great success recent years. However, there is still a 
vast interest beyond standard model. It is well known that the process b — > sy is extremely 
sensitive to new physics beyond the Standard Model||. It has been argued that this experiment 
provides more information about restrictions on the Standard Model, 2-Higgs doublet model, 
Supersymmetry, Technicolor and etc. Since the top quark is much heavier than other fermions, 
and its interactions may be quite sensitive to new physics, the interactions of the top quark 
are of special importance. In ref.0, a right-handed coupling of the top-bottom charged current 
is added to the standard model(SM). Based on this, the authors give out the constraints to 
right-handed top quark current by b — > 57 decay. In fact, the right handed current is also a low 
energy phenomena of SU(2)i x SU(2)r x U(l) left-right symmetric model||. In this model, 
the right handed charged current is induced by the Wl — Wr mixing in addition to the usual 
left-handed current in the standard model. Interference of the right- and left-handed currents 
in the penguin diagram enables a chirality-flip by the top quark mass m t inside the loop, and 
therefore leads to an amplitude proportional to m t || [|, |J . (In standard model, the amplitude is 
proportional to the mass of the bottom or strange quark, because the pure (V-A) structure of the 
charged currents requires the chirality-flip to proceed only through the mass of the initial or the 
final state quark.) Therefore, large contributions from top-bottom quark right-handed current 
to b — > 57 amplitude occur. To first order approximation, this just like to add a right-handed 
current to top-bottom quark coupling in SM. 

Recently the CLEO Collaboration has measured the inclusive branching ratio of b — > 37 to 

be0 

Br{b -> sy) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10~ 4 . (1) 

Corresponding to 95% confidence level, the range is 1 x 10" 4 < Br(b -y 57) < 4 x 10~ 4 . 
This is a more stringent constraints compared with previous CLEO experiments || cited by the 
above papers [|], f|, (5|. Furthermore, with more precise experiments, a more accurate theoretical 
calculation of this decay rate is also needed. 



2 



The inclusive b — > s'y decay rate is often assumed to be well described by the spectator 
model, where the b quark goes a radiative decay The QCD corrections to this decay have been 
calculated to leading logarithmic accuracy by many authors || [TT| O, [TJ, [T3], FJJ], and are known 
to enhance the decay rate within SM by a factor of 3-4. This enhancement, however, is subject 



to large uncertainties due to many reasons One of these reasons is the QCD running from 
top quark scale to W boson scale, which has been discussed in SM by ref.|L6|, O, in 2-Higgs 



doublet model by ref.JT^], and in Supersymmetry case by ref.|19[. This contribution which is 
usually considered as a next-leading order effect was found to give an additional enhancement 
up to 20% in the case of a much heavy top quark mass. 

In our present paper, by introducing a right-handed charged current to top-bottom couplings, 
we recalculate the b — > 37 decay including QCD running from m top to M w , in addition to 
corrections from My/ to m^, in order to give a more complete leading log result in this model. 
Since the branching ratio is proportional to |/jfm t /mb| 2 , it is more complicated than that in SM 
or other models. Finally more recent CLEO experiment is used to give more precise restrictions 
to this anomalous top quark coupling. 

In the next section, we first integrate out the top quark, generating an effective five-quark 
theory. By using the renormalization group equation, we run the effective field theory down to 
the W-scale where the weak bosons are removed. Then we continue running the effective field 
theory down to b-quark scale to include QCD corrections from Mw to m^. In section 3, the rate 
of radiative b decay is obtained. Section 4 is a short summary. 



2 QCD Corrections to b —> 57 decay 

In standard model and many other models, there is no flavor changing neutral current at tree 
level. It is only occurred through electroweak loop. In standard model Lagrangian, the relevant 
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charged current in gauge reads 



s 



Wt 



+ 



1 u e/2 



V 92 



V2 M w 
+ h.c. 



u c t ) M V V 

R 



S 

\ h J 



u c t )VM D 



s 

\ b I 



R- 



(2) 



Where V represents the 3x3 unitary Kobayashi-Maskawa matrix, Mjj and Mp denote the 
diagonalized quark mass matrices, the subscript L and R denote left-handed and right-handed 
quarks, respectively. This charged current, which is the main cause of flavor changing neutral 
current, is pure left-handed. 

Like ref. |IJ, a phenomenological coupling of the right-handed top and bottom quarks to the 
W boson is introduced: 



£ = g2/V2^ 2 V tb f%t Rl ,b R W£ + 



V tb f R b t(m t P R - m b P L )b<P + + h.c. 



(3) 



V2M W 

with f R b denoting the strength of this additional coupling. Using this interaction (^) together 
with (H), we performed the QCD corrections to b — >• 37 decay. 

At first, we integrate out the top quark, generating an effective five quark theory, intro- 
ducing dimension-5 and dimension-6 effective operators as to include effects of the absent top 
quark. Higher dimension operators are suppressed by factor of p 2 /m 2 , here p 2 characterizes the 



interesting external momentum of b quark p 2 



~ mi 



For leading order of ml/m 2 , dimension-6 



operators are good enough to make a complete basis of operators |L6], |T7| 3 : 



dimension 5 : 
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; m b s L D 2 b R , 



3 Notice here W^ R is a new operator compared with SM case. 



0\ R = ^ 2 ^- 2 m b s L a^X% R G% 

0\ R = ^ 2 ^m b s L a^b R F^ 

Qlr = H € glm b (f) + (f)_s L b R , 

W l LR = -t^g 2 3 m b WlW^s L a^b R , 

Wlr = V e gim b W!£W?s L b R , 
dimension 6 : 

Pl A = -^s L T* a D»D»D°b L , 

Pi = I^^l^F^ 

Pi = i^^Ks^D^ 
lbir z 

R 1 l = in e gi<j)+<j)-s L pb L , 
R 2 L = ifj, e gl(D a (j) + )(j)Js L 'y a b L , 

W\ = iii e glWlW^ Ll ,ip lv b L , 
W 2 L = %^gl[D^WX)W^ Llill(ylv b Ll 
Wl = zSg 2 3 W + ,W»s L flb L , 

W A L = w e glW»W»s L (D^ v + ^D„)b L . (4) 
Where siD^lubi stands for {shD^ubL + (D^s^^bi) and the covariant derivative is defined as 

D„ = d„- z^ 2 g 3 X a G a u - i^ 2 eQA„ 



with g 3 denoting the QCD coupling constant. The tensor T A appearing in P^ A assumes the 
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following Lorentz structure, the index A ranging from 1 to 4: 

(5) 



Pjj.ua ~ 9valjli Pjj.ua ~ ^jivarl 75- 

In our following calculations, we try to work in a background field gauge (with £ 



1)[ 20 1, in order to maintain explicit gauge invariance in calculations of anomalous dimensions. 



Furthermore the usually trilinear interaction between photon, W boson, and would-be Goldstone 



boson vanish in this gauge. So that, we did not include operators which involve both W boson 
and would-be Goldstone boson in one operator, like W+<\)-~s-i,^ v b L . With the above operators, we 
can write down our effective Hamiltonian as 

H eff = 2V2G F V tb V t l £ (6) 

i 

The coefficients Ci(m t ) of operators can be calculated from full theory by matching conditions, 



keeping only leading orders of p /mffllE, I7| . Terms proportional to f^mt/rrib are from right 



handed current. Since mt/m& is a very large value, it is very convenient for us to keep only 
leading orders of m t jm b . This is also the reason why only top-bottom right-handed charge 
current is introduced, while other up-down, charm-strange right-handed currents are ignored. 
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(7) 

The other coefficients of operators are all from the integrations of electroweak loops. Terms like 
log(/i 2 /m 2 ) in coefficients of operators 0\ R , P| and Pi vanish here, because of the matching 
scale fi = m t . They will be regenerated at lower scales by renormalization group running of 
electroweak later. The other logarithms are all from the finite part integration of loops, for there 
are two different mass scale particles in one loop. 

f 1 + 5 5 . _\ £tb m t fl-35-45 2 5 - A5\ \ 
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where 5 = M^/m 2 . When /f = 0, the above results (0), (g) reduce to that of SM case[[L7|. 
The renormalization group equation satisfied by the coefficient functions Ci(fi) is 

^ C *V) = E(7 T )^»- 



(9) 



The solution to this renormalization group equation (||) appears in obvious matrix notation as 
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(10) 



where the anomalous dimension matrix 7^ is calculated in practice by requiring renormalization 
group equations for Green functions with insertions of composite operators to be satisfied order 
by order in perturbation theory [|16| , |T7j . 
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These mixings are all between operators (Q,R,W) induced by tree-diagram and operators (0,P) 
induced by loop-diagram. In order to see how the renormalization group method is accomplished, 
we neglect the proper QCD corrections for the moment, so that we can take into account only 



the above entries of anomalous dimensions. Insert this matrix to eqn.(|10|), we find the following 
relations: 



c o1r {m w ) 

C Pl (M w ) 
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C ol (mj + fi—f&log^, 
LR nib irif 

C P 2 (mt) - \ log - 5 log , 
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Here the renormalization group equation reproduces the log(/i 2 /mf) terms in the coefficients 
of operators at equation (H) which vanished at /x = m t . This proves the consistence of the whole 
calculation. 

The QCD anomalous dimensions for each of the operators in our basis are [[16], 17 1: 
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with such complicated anomalous dimensions with so many operators, the algebra computation 
code REDUCE is used. Combining the above three matrices ( |TT| , |T2],[I3]) into one large 18x18 



matrix, we diagonalize it analytically to give 18 eigenvalues and 18 eigenvectors. Using equation 



fliTf ), one can have the coefficients of operators at \i = M\y- All analytic expressions of coefficients 
at this scale, which are functions of coefficients at fi = m t , are included in the appendix. 

In order to continue running the basis operator coefficients down to lower scales, one must 
integrate out the weak gauge bosons and would-be Goldstone bosons at /i = Myy scale. From 
matching conditions, one finds the following relations between coefficient functions just below(-) 



and above(+) /i = Miy [fLo] , |T7jj : 
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+ 2/9, 



C p i4Mw) = C p i, 4 (M+) + l, 
C Pl (M w ) = C Pl {M+)-C wl {M+)- 3/2, 
C Pl (Mw) = C P 3(M+), 

C Pt (Mw) = C P 4(M+) + 9. (14) 
The operators involving Goldstone cf) and W bosons are absent after this matching. 

In addition to these, there are new four-quark operators from the matching|9|, 12], |14| : 

q q 

with coefficients 

Ci(M w ) = 0, 2 = 1,3,4,5,6, C 2 (M W ) = 1. 

Although only coefficient of 2 is not zero, the other operators should also be included, for their 
mixing through renormalization group running. 

To simplify the calculation and compare with the previous results, equations of motion(EOM) |21] 
is used to reduce all the remaining two-quark operators to the gluon and photon magnetic mo- 
ment operators 0\ R and 0\ R . Neglecting the strange quark mass in comparison with the bottom 
quark mass, we obtain the on-shell equivalence relations: 

O 1 = Pr 1 ' 1 = - -O 2 - -O 3 

LR L ^ LR ^ LR) 

Pl = ~Pl = ~ — 4^^' 

Pi' 3 = P 2 L = 0, 

Pi = -\01r- (16) 



To be comparable with previous results without QCD corrections from m top to Myy, opera- 
tors 3 LR , 0\ R are rewritten as 7 , 8 like ref.[|, |TJ, [TJ]], 

7 = (e/16n 2 )m b s L cr^b R F^, 

8 = (g/167T 2 )m b s L a^T a b R G" (17) 
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Then the operator basis now consists of 8 operators. Using eqn.flliD ( Ti6|) together with ( p4] ) in 
the appendix, the effective Hamiltonian appearing just below the W-scale is easily drawn out: 
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a s (M w ) 

Since they are expressed by coefficients of operators at fi = m t and QCD coupling a s , it is 
convenient to utilize these formulas. 

If QCD correction from m t to M w is neglected[by setting a s (M w ) = a s {m t ) in the above 
equations (0),(^3)], the results are reduced exactly to those the previous authors gave|], [|. 

The QCD corrections from \i = M\y to fi = nib have been a subject of many papers|| 



10], [11], [12], [13], [14]]. It has attracted great interest recent years, since most QCD corrections 
come from this stage. For a long time, there are some discrepancies in this calculation, for it is 
rather lengthy calculation involving 2-loop diagrams. Until recently, it is completely resolved, 
and a complete leading log result is given fl4]j. Here, we are free to utilize these results, so the 
coefficients of operators at /i = nib scale are 



c 7 e// K) = v 16/23 c 7 (m w ) + - ^)C,{M W ) + C 2 (M W ) £ h*p 

6 i=l 



(21) 
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With r] = a s (M w )/a s (m b ), 



/626126 56281 3 1 > 
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V272277 51730 7 14' ' ' j 
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Here m\> = 4.9GeV is used. 



3 The B —> X s 7 decay rate 

Following previous authors^, |12|, [14||, one obtains the B — > decay rate normalized to the 
quite well established semileptonic decay rate BR(B — > X c eu). 

BR(B -> X s -f)/BR(B -> X c eV) ~ r(6 -> s-y)/T(b -> cez7). (22) 

Then 

- -^^ICJ'W f 1 - ^/KW) , (23) 
BR(B —> X c eis) TTg(m c /m b ) \ Sn J 

where g(m c /mb) — 0.45 and f(m c /mb) — 2.4 correspond to the phase space factor and the 
one-loop QCD correction to the semileptonic decay, respectively^]. The electromagnetic fine 
structure constant evaluated at the b quark scale takes value as ctQED{ m b) — 1/132.7. If we take 
experimental result BriB — > X c eV) = 10.8% [Q, the branching ratios of B — > X^ is found. 
Compared with the previous results without QCD running from m t to M^[0, the branching 
ratio is more sensitive to the value of It is no longer a pure enhancement like SM case. It is 
rather suppressed when \ f^\ takes large value, for example, When \ fj^\ = —0.08, the branching 
ratio is only about 1/10 of that with QCD correction from m t to My/ neglected. While when \f^\ 
takes small value, the branching ratio is rather enhanced, especially when |/*| = 0, it is enhanced 
11% corresponding to SM case. Since the right-handed current contribution is proportional to 
|/Jj b | 2 , the whole QCD correction to b — ► s'j decay is always a large enhancement. 

In Fig 1, the branching ratio of b — > s'j is displayed as functions of f§ with mt = 174GeVp4|. 



The three lines correspond to a s {mz) = 0.107, 0.117, 0.127. They are all parabolic lines, because 
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Br{b — > 57) is proportional to |CV(m&)| 2 , and C^im^) is proportional to /Jj. That the three lines 
are very close to each other implies that the branching ratios are not sensitive to the values of 
a s . 

Fig.2 is branching ratios displayed as functions of /jf, with different top quark masses 
158GeV, 174GeV and 190GeV. The QCD coupling constant is a s {m z ) = 0.117. Since contribu- 
tions from right-handed current are proportional to \f^m t /mb\ 2 , the three lines in Fig.2 diverge 
in the region far from = 0, but nearly converge to a same line when /Jf — > 0. 

With the recent CLEO experiments of b — > 37 decay, a new constraint to right-handed 
current can be found. From Fig.l and Fig.2, the parameter is constrained in two small 
windows -0.087 < /f < -0.050 and -0.023 < /f < 0.002. Compared with results without 
QCD corrections from mt to Mw0, the whole line transfers to left side(towards minus). So the 
windows also transfer to minus side. Most allowed values of /Jf are negative. 

4 Conclusion 

In the above, we have introduced an anomalous right-handed current to top, bottom quarks and 
W boson coupling in the SM, and given the full leading log QCD corrections (including QCD 
running from m top to M w ) to b — > 37 decay in this model. 

The QCD corrections from m top to M\y to b —>■ 57 decay enhance the decay rate in small 
values of like in SM case; but suppress the decay rate when takes larger value. The 
whole QCD correction still makes a large enhancement. The restrictions from b — > 57 decay to 
anomalous top quark coupling are strict, only two narrow windows are allowed. It is shown that 
the decay b — > S7 has been the most restrictive process so far in constraining the parameters of 
the right-handed current of top quark. 

If a complete QCD next-leading log result of b — >• S7 decay is performed, one can expect to 
obtain more precise results from b — > 37 decay by freely using our above results. 
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A Operator coefficients at fi = scale 

We quote here the results of the coefficients of operators at fi = M^> scale. The W bosons and 
would-be Goldstone bosons are not integrated out yet. 

c o1r (m w ) = (2C 14 / 23 -C 8/23 )^K) + 4(-C 14/23 + C 8/23 )^o L K) 
+ (2C 14/23 -C 8/23 -C 4/23 )^mK) 

f c 14/23 _ ^9_ C 8/23 _ US 4 /23 864 Q ( ) 

V 130 s 274 s 8905 s / p l v * ; 

A_ C 14/23 + Hi C 8/23 _ 113 C -4/23 + ^1^13/138^) Q ( } 

V 130 s 274 s 8905 s / p l V ; 
= \ (C 14 / 23 - C 8/23 ) C ol R (m t ) + (-C 14 / 23 + 2C 8 / 23 ) C olR (m t ) 

+^(( 14/23 -(- 4/23 )C p uW 

+ fI C 14/23 _ A C 8/23 _ 113^-4/23 + i 3 ^ 1 13/138 \ 

V4 65 s 548 s 8905 s / p ^ V ; 

+ (_I C 14/23 + i3_ C 8/23 _ 113 c _ 4 /23 + 1 13/1! «A ^ ( } 

V 4 s 130 s 548 s 8905 s / p l V ; 
C s fl (M w ) = (5C 14/23 - ^C 8/23 - ^C 16/23 ) C olR (m t ) + (-10C 14/23 + 2C 8 / 23 + 8C 16 / 23 ) Q 

+(^ 3 C olR (m t ) + (sc 14 / 23 - \c A/2Z - ^C 16/23 ) C Pl 4m) 

, (\ .14/23 _ _^_a8/23 _ I^a-4/23 , 43 2 ^n 3 /i 38 9 l 6/23 \ 

+ U C 65 C 548 C + 8905 C 4 C J G ^ 2 ^ 

^5,14/23 + j^.8/23 _ 113^-4/23 + ^,113/138 + ^16/23^ C Am t ) 

V 2 s 130 s 548 s 8905 s 4 s / p ^ V ; 
+\ (l " C 16/23 ) C Pt {m t ) + g (i _ C -i) 167r 2 C^(M w ), 

c p m(m w ) = c P ,a = r 4/23 c P MK) + §^ (C 4/23 - C 113/138 ) (CpuW + C p mH), 
^m(m^) = (C 8/23 - C 4/23 ) C Pl 4m) + Qc 8/23 - ^jc 4/23 + ^|C 113/138 ) C Pl ,(m t ) 

-\( 8/23 - ^C 4/23 + ^|C 113/138 ) C P M(m,), 



+ - 
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c Pl 4M w ) = (-C 8/23 + C 4/23 ) c P uH + (-ic 8/23 + ^jc 4/23 + ^C 113/138 ) c Pl 4m) 
+ Qc 8/23 + ^C" 4/23 + ^C 113/138 ) C Pl 4m), 

C Pl (M w ) = ^(c 113 / 138 -l)(C p i, 2 (m t ) + C p MK))+C P| K) 



Here ( = a s (m t )/a s (M w ). The factor of IQir 2 in some of these expressions arises from mixing 
between operators induced by tree diagrams and those by loop diagrams. 
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Figure Captions 

Fig.l BR(B -> X s7 ) as function of /f for different a s values. a s {m z ) = 0.107, 0.117, 0.127. 
Fig. 2 BR(£> — > X s 7) as function of /jj for different top quark masses. m t = 158, 174, 190GeV. 
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